Quantum error correction of coherent errors by randomization 
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A general error correction method is presented which is capable of correcting coherent errors 
originating from static residual inter-qubit couplings in a quantum computer. It is based on a 
randomization of static imperfections in a many-qubit system by the repeated application of Pauli 
operators which change the computational basis. This Pauli-Random-Error-Correction (PAREC)- 
method eliminates coherent errors produced by static imperfections and increases significantly the 
maximum time over which realistic quantum computations can be performed reliably. Furthermore, 
it does not require redundancy so that all physical qubits involved can be used for logical purposes. 
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Current developments in quantum physics demon- 
strate in an impressive way its technological potential 
0. In quantum computation, e.g., characteristic quan- 
tum phenomena, such as interference and entanglement, 
are exploited for solving computational tasks more ef- 
ficiently than by classical means 0, 0, 0, IE El ■ How- 
ever, these quantum phenomena are affected easily by 
unknown residual inter-qubit couplings or by interactions 
with an uncontrolled environment [7|. In order to pro- 
tect quantum algorithms against such undesired influ- 
ences powerful methods of error correction have been de- 
veloped over the last years. 

So far techniques of quantum error correction have 
concentrated predominantly on decoherence caused by 
uncontrolled couplings to environments 0, 13 ■ I* 1 t nese 
cases appropriate syndrome measurements and recovery 
operations can reverse errors. However, up to now much 
less is known about the correction of coherent, unitary 
errors. Even if a quantum information processor (QIP) is 
isolated entirely from its environment and if all quantum 
gates are performed perfectly, there may still be resid- 
ual inter-qubit couplings affecting its performance. Re- 
cently, it was demonstrated that static imperfections, i.e. 
random inter-qubit couplings which remain unchanged 
during a quantum computation, restrict the computa- 
tional capabilities of a many-qubit QIP significantly as 
they cause quantum chaos and quantum phase transi- 
tions |ld| . Furthermore, in addition to a usual exponen- 
tial decay such static imperfections also cause a Gaussian 
decrease of the fidelity with time. At sufficiently long 
times this Gaussian decrease dominates the decay of the 
fidelity thus limiting significantly the maximum reliable 
computation times of many-qubit QIPs |IU Il2| . 

In this Letter a general error correcting method is pre- 
sented for overcoming these disastrous consequences of 
static imperfections. It is based on the repeated random 
application of Pauli operators to all the qubits of a QIP. 
The resulting random changes of the computational ba- 
sis together with appropriate compensating changes of 
the quantum gates slow down the rapid Gaussian decay 
of the fidelity and change it to a linear-in-time exponen- 



tial one. As a result this Pauli-Random-Error-Correction 
(PAREC)-method increases significantly the maximum 
time scale of reliable quantum computation. In addition, 
neither control measurements nor redundant qubits are 
required so that all physical qubits are logical qubits. 

In order to put the problem into perspective let us 
concentrate on the quantum algorithm of the quantum 
tent-map as a particular example [p3 |. One iteration of 
this special case of a quantum-rotator-map is governed 
by the unitary operator 



U = 



-iTp 2 /(2mh) -ikV(x)/h 



(1) 



It describes the one-dimensional dynamics of a periodi- 
cally kicked particle of mass m. The operators p and x 
denote momentum and position operators and T is the 
period of the kicks of magnitude kV(x). The dynamics 
of the particle is assumed to be confined to the spatial 
interval < x < I with periodic boundary conditions. 
The name of this quantum-map originates from the force 
which resembles the form of a tent, i.e. 



V'(x) 
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(0 < x < 1/2) 
(1/2 < x < I). 



(2) 



Due to the periodic boundary conditions the momentum 
eigenvalues are given by p n = 2-Khn/l with n £ Z. Impos- 
ing the 'resonance condition' T = [m(l / '2tt) 2 / 'H\(2ir / 'N) 
with N £ N implies the symmetry (p n +N I U \ Pu'+n) = 
(—l) N (p n | U | p n ') and U decomposes into a direct 
sum of N x N matrices 13]. Thus, for a given value of 
N the dynamics of the quantum tent-map can be simu- 
lated on a quantum computer (QC) with n q qubits pro- 
vided N = 2 nq . In this case the unitary operation of 
Eq.Q can be performed with the help of n g — (9/2)n^ — 
(\l/2)n q + 4 universal quantum gates, i.e. Hadamard-, 
phase-, controlled-phase-, and controlled-not gates [l^ . 
The classical limit of the quantum tent-map corresponds 
toT^ 0,/c-> oo with K = kT[l/(2TTH)]/[m(l/2TT) 2 /h] 
remaining constant. In this parameter regime the tent- 
map exhibits all complex dynamical features character- 
istic for quantum chaos . 
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In order to model static imperfections we assume that 
the n q qubits of a realistic QC are coupled by ran- 
dom Heisenberg-type nearest-neighbour interactions as 
described by the Hamiltonian 

n q —l n q —2 

1 + YiYi+l + ZiZi + i) 

t=0 i=0 

(3) 

with the Pauli (spin) operators X,Y,Z. The quan- 
tities Si and Ji denote the strengths of the detuning 
and of the nearest-neighbour interaction of qubit i. In 
the case of static imperfections these quantities are dis- 
tributed randomly and homogeneously in the energy- 
interval [— \firi, yf?>rj\ and remain static in time during 
a quantum computation. 

In general, after t iterations of a quantum map the fi- 
delity, defined through the ideal and the perturbed quan- 
tum states | ip(t)) and | tl) n (t)), i.e. / =| (il>(t) | ip v {t)) 
decays according to [l2] 



ln/(t) 



tr. 



t c tH 



(4) 



This relation is valid as long as rj and t are sufficiently 
small so that the fidelity / remains close to unity. The 
time scale t c governing the linear-in-time exponential de- 
cay is determined by Fermi's Golden rule. In partic- 
ular, in recent simulations 12] it was found to be in- 
versely proportional to (n q n^rj 2 ). The second character- 
istic time scale entering Eq.Q is the Heisenberg time 
t H w 2"" which is determined by the dimension of the 
n g -qubit Hubert space. According to Eq.Q the linear- in- 
time exponential decay changes to the much more rapid 
quadratic Gaussian decay roughly after t sa t# iterations. 
Such quadratic-in-time Gaussian decays are characteris- 
tic for coherent dephasing phenomena. The quadratic 
fidelity decrease of Eq. (0} is drastically limiting the max- 
imum time over which a quantum computation can be 
performed reliably. Contrary to static imperfections, ran- 
dom imperfections which change from gate to gate, e.g., 
lead to a purely linear-in-time exponential decay of the 
fidelity [ll|- This suggests the idea that a randomiza- 
tion of static imperfections might help to slow down the 
fidelity decay to a linear-in-time one thus increasing sig- 
nificantly the reliable computation time. 

But how can a randomization of static imperfections be 
achieved efficiently? A basic idea of quantum error cor- 
rection is to exploit the freedom of choice of the compu- 
tational basis for an appropriate encoding. This idea can 
also be used for an efficient randomization of static im- 
perfections by changing the computational basis repeat- 
edly and randomly during a quantum computation. How- 
ever, to leave the quantum algorithm unchanged these 
basis changes have to be compensated by appropriate 
transformations of the universal quantum gates. 

In order to address this issue let us concentrate on 
the particular example of the quantum tent-map algo- 
rithm and on the static Heisenberg-type imperfections 
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FIG. 1: The basic idea of the PAREC-method: The two boxes 
(full lines) represent two sequences of universal quantum gates 
for n q = 4 qubits. Two random sequences of Pauli operators 
(Xi, Y2, Z3, I4) and (Y, X2, X3, Z4) are also indicated. The 
unitary Pauli operators outside the dashed boxes (full lines) 
are applied to the qubits whereas the ones inside the dashed 
boxes (dashed lines) are taken into account by appropriate 
permutations of the elementary quantum gates. Due to the 
identities X 2 = Y 2 — Z 2 — 1 the inserted random sequences 
of Pauli operators change the computational basis but leave 
the ideal quantum algorithm unchanged. 



described by Eq.JJJl. A convenient way of realizing such 
random changes of the computational basis is to apply 
repeatedly randomly selected Pauli operators to all the 
n q qubits of the QC. For this purpose it is advantageous 
to represent the elementary quantum gates of the quan- 
tum algorithm in terms of a special Hamiltonian set of 
universal quantum gates, namely 



S ±X] (A<p) = 
S ±XkX] {A<f) 



S±z 3 (Acp) 



(5) 



The Hamiltonians appearing in the exponents of Eq.JSJ) 
are themselves Pauli operators. Therefore, any unitary 
transformation R originating from Pauli operators ei- 
ther leaves these Hamiltonians invariant or changes their 
signs, such as 



s ±Xj (A^) if ki e {i 3 ,x 3 } 
S TXj (A<t>) if RjeiY^Zj}. 



(6) 



Thus, with the help of these Hamiltonian quantum gates 
any change of the computational basis originating from 
a randomly selected set of Pauli operators can be com- 
pensated by an appropriate permutation of the universal 
quantum gates of Eq.©. 

On the basis of these considerations during a quantum 
computation an efficient correction of static errors can 
be achieved by the PAREC-method in the following way 
(compare with Fig.l). In the first step randomly selected 
unitary operations from the set {1, X, Y , are applied 
to all n q qubits of the QC, say (X%, Y2, )• The infor- 

mation about which qubit has been transformed by which 
Pauli operator is stored in a classical memory. In the sec- 
ond step one starts the quantum computation by apply- 
ing a sequence of properly permuted universal quantum 
gates (compare with Eqs.© and iJBJ) as described by the 
first dashed box of Fig.l. This simple permutation does 
not require any extra significant computational effort. In 
the third step a second sequence of Pauli operators is se- 
lected randomly, say (Zi, Y~2, ...,X nq ), and the combined 
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quantum gates (XiZi, Y2Y2, —, In q ^n q ) are applied to 
all the n q qubits of the QC. These combined quantum 
gates are again Pauli operators. The information about 
the spin operators of the second selection is again stored 
in a classical memory. Afterwards the second sequence 
of properly permuted universal quantum gates is per- 
formed (second dashed box of Fig.l). In the subsequent 
stages of the PAREC-method these steps are repeated 
after sequences of universal quantum gates of appropri- 
ate lengths n go f . The influence of the choice of n ge f on 
the error correction will be discussed later (compare with 
Fig. 3). Finally, after the application of the last quantum 
gate the last randomly selected sequence of Pauli opera- 
tors is applied to all qubits. As apparent from Fig.l this 
PAREC-method leaves the ideal quantum algorithm un- 
changed. However, the repeatedly applied random uni- 
tary transformations produced by the Pauli operators 
change the signs of the parameters Si and Ji of Eq.© 
thus causing a randomization of the static errors. As a 
result we expect a significant improvement of the fidelity 
decay. 
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FIG. 2: (Color) Quantum Poincare sections with Husimi- 
functions at t — 3000 in scaled momentum and position vari- 
ables y =p[(l/(2nh)] G [0, 2vr] and x = x[2tt/1] £ [0,2tt]: The 
parameters are K = 1.7 and n q = 10. The initially prepared 
coherent states are centered around (tt/4, 0) (left panel) and 
(5.35, 0) (right panel). First row: ideal dynamics; second row: 
static imperfections with e = 5 x 10~ 6 ; third row: PAREC- 
method applied after each sequence of n gc f = 20 universal 
quantum gates of Ref. |T^|. The probability density is coded 
in colors (red/maximum, blue/zero). 

The stabilizing properties of the PAREC-method are 




t/1000 ln(t) 



FIG. 3: Dependence of the fidelity f(t) on the number of it- 
erations: Parameters as in the left panel of Fig. 2; left: static 
imperfections without error correction, PAREC after each 
map iteration, after each n ge f = 50, and after each n gc f = 20 
quantum gates (from bottom up); right: static imperfections 
without error correction, with PAREC after each map iter- 
ation and after each sequence of n gc f = 20 quantum gates 
(full curves), best fits for linear- and quadratic-in-time decays 
(dashed curves). 

investigated numerically in Figs. 2 and 3 where it is ap- 
plied to the iterative quantum tent-map. For this pur- 
pose the unitary operator of Eq.Q is decomposed into 
a sequence of h g universal quantum gates of the form of 
Eq.(J5}. This can be achieved in a straightforward way, 
e.g., by corresponding replacements of the already known 
gate decomposition into n g gates [l^- The PAREC- 
randomization is applied after appropriately chosen se- 
quences of n ge f quantum gates of this latter decomposi- 
tion of n g gates. The influence of static imperfections 
is modeled by assuming that these universal quantum 
gates are performed instantaneously but that there is 
a certain time delay between any two successive quan- 
tum gates during which static imperfections cause er- 
rors. This time delay At models in an approximate way 
the time required for readjustments of the control unit 
of the QC before it can control the next quantum gate. 
As the quantum gates of Eq.JSJ) involve different accu- 
mulated phases Acj>, within our model we also assume 
that quantum gates with larger phases require longer 
readjustment-times. Correspondingly, the readjustment- 
times At\ and At2 after two successive quantum gates 
with accumulated phases A<pi and Afa are related by 
At 2 /Ati — A<f>2 / A<f>\. Thus, for a typical sequence of 
two quantum gates, e.g., the influence of static imperfec- 
tions is modeled by 

... [e -i(HAt/H)(At n+1 M SxiXj{A(l)n+i)] 

[ e -HAAt/n)(MnMs Xk ^)] ... (7) 

with At denoting the readjustment-time associated with 
a phase change | A<fi |= tt. Correspondingly, the param- 
eters characterizing the average strength of the influence 
of the static imperfections between successive quantum 
gates are (SiAt/h) and (JiAt/K) which are selected ran- 
domly and uniformly from the interval [— v3e, V3e] and 
which remain constant during the quantum computation. 

The ideal dynamics of the quantum tent-map is illus- 
trated in the first row of Fig. 2 where its Husimi-functions 
[l3| are plotted after t = 3000 iterations for two initially 
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prepared coherent states located close to the classically 
unstable (left) and to the classically stable (right) fixed 
points. In the left figure of the first row the initially 
prepared coherent state spreads almost uniformly over 
the classically chaotic component of phase space. Clas- 
sically inaccessible regions originating from Kolmogorov- 
Arnold-Moser-tori are also apparent. The quantum prob- 
ability leaking into these regions by quantum tunneling 
is still negligibly small. In the right figure of the first row 
the quantum probability is still concentrated in the re- 
gion of the initially prepared coherent state. This reflects 
the approximately regular classical dynamics in this part 
of phase space. Static imperfections modify these dy- 
namical characteristics significantly, as is apparent from 
the second row of Fig. 2. Most prominently the influence 
of quantum tunneling into the classically inaccessible re- 
gions of phase space is no longer negligibly small. Fur- 
thermore, the detailed structures in both the regular and 
the chaotic parts of phase space are modified significantly. 
The corresponding results of the PAREC-method are de- 
picted in the third row. Despite the fact that random 
sequences of Pauli operators are applied only after each 
sequence of n gc f = 20 gates of Ref. 0] the corrected 
quantum states resemble the ideal case very closely. 

The quantitative dependence of the fidelity on the 
number of iterations of the quantum tent-map is de- 
picted in Fig. 3 for the same parameters as in the left 
panel of Fig. 2. The initially prepared coherent state of 
the right panel of Fig. 2 yields similar results. It is appar- 
ent that the PAREC-method changes the fidelity decay 
from quadratic- to linear-in-time, i.e. In /(f) = —t/t c . 
Thus, the time over which quantum computations can 
be performed reliably is increased significantly. In ad- 
dition, one also notices that with increasing numbers of 
random unitary operations the decay rate l/t c also de- 



creases. In particular, the best fits to the PAREC-results 
of Fig. 3 together with further numerical studies suggest 
a dependence of the form l/t c = ae 2 n q n g ng e [ as long as 
1 < n gcf ^ n g- Thereby, Tigcf is the effective number 
of original gates of Ref. [12j over which the influence of 
static imperfections adds up coherently. Our numerical 
data of Fig. 3 give a ~ 1. Physically speaking this depen- 
dence is plausible on the basis of the following heuristic 
consideration. If the PAREC-method is repeated after 
each map iteration, one observes n ge f = n g as coherence 
is destroyed by the random basis changes and by the 
chaotic dynamics after each iteration. This result is con- 
sistent with numerical studies |12| . In the extreme oppo- 
site case in which the PAREC-method is repeated already 
after each universal quantum gate one expects n gc f = 1 
as coherence is destroyed by random basis changes al- 
ready after each gate operation. The above mentioned 
dependence interpolates linearly between these two ex- 
treme cases. However, due to fluctuations it is expected 
that these considerations only apply for sufficiently large 
values of n go f. 

In summary, a general method for the correction of uni- 
tary static inter-qubit errors has been presented. This 
PAREC-method is particularly well suited to stabilize 
many-qubit systems against the disastrous effects of 
static imperfections in arbitrary quantum algorithms. In 
contrast to conventional quantum error correcting meth- 
ods which exploit redundancy the PAREC-method does 
not require any extra qubits so that all physical qubits 
can be used in an optimal way. 
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